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ABSTRACT 

Uncertainty in the estimate of the microwave propagation delay through the troposphere is a 
principal limiting factor to the accuracy of the technique of very long baseline interferometry 
(VLBI). This uncertainty is due primarily to tropospheric water vapor, the total amount and verti- 
cal distribution of which is variable. Because water vapor both delays and attenuates microwave 
signals, the propagation delay, or wet path length, can be estimated from the microwave brightness 
temperature near the 22.235 GHz transition of water vapor. 

We analyzed the data from a total of 240 radiosonde launches taken simultaneously in 1974 at 
Chatham, Massachusetts; Albany, New York; and Portland, Maine. Estimates of brightness tempera- 
ture at 19 and 22 GHz and wet path length were made from these data. The wet path length in the 
zenith direction could be estimated from the surface water vapor density to an accuracy of 5 cm for 
the summer data and 2 cm for winter data. Using the brightness temperatures, the wet path could 
be estimated to an accuracy of 0.3 cm. 

Two dual-frequency radiometers constructed by the National Radio Astronomy Observatory 
(NRAO) were refurbished in order to test these techniques. These radiometers were capable of 
measuring the difference in the brightness temperature at 30° elevation angle and at the zenith to an 
accuracy of about 1°K. In August 1975, 45 radiosondes were launched from Haystack Observatory 
over an 1 1-day period. Brightness temperature measurements were made simultaneously at 19 and 
22 GHz with the NRAO radiometers. The rms error for the estimation of wet path length from 
surface meteorological parameters was 3.2cm, and from the radiometer brightness temperatures, 
1.5 cm. 
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INTRODUCTION 

The neutral atmosphere retards and attenuates propagating electromagnetic waves. The retardation 
limits the accuracy to which very long baseline interferometry can be used to measure radio source 
positions and baseline vectors between the antennas. Atmospheric water vapor can contribute up to 
about 40 cm of excess propagation path length in the zenith direction at microwave frequencies. 
The exact amount cannot be predicted accurately from ground level meteorological variables since 
the water vapor is not well mixed in the atmosphere. However, since the index of refraction and ab- 
sorption coefficient are functions of the water vapor density, the brightness temperature due to the 
self-emission, a weighted integral of the absorption coefficient, is related to the path length, the 
integral of the index of refraction. Hence, ground-based radiometric measurements can be used to 
estimate the excess phase path. Early evaluations of the effectiveness of this technique were made 
by Waters (1967) and Shaper, Staelin, and Waters (1970). 

The vertical profiles of temperature and water vapor density are routinely obtained from radio- 
sondes launched by the National Weather Service (NWS). From these data, the index of refraction 
can be calculated and the path delay estimated. The microwave absorption coefficient can also be 
calculated and the brightness temperature estimated from the equation of radiative transfer. We 
completed a theoretical study based on radiosondes launched at three locations in New England in 
order to determine how well the path delay could be estimated from: (1) surface meteorological 
data, (2) radiosonde data from a remote station, and (3) microwave radiometry. We also conducted 
an 1 1-day experiment during which we compared actual radiometry data from two microwave radi- 
ometers with values of brightness temperatures and path length derived from radiosonde data. Our 
study was focused on establishing how well various techniques work. We used the method of linear 
regression analysis and examined the residuals. Substantial progress has been made recently by Wu 
(1979) and Claflin, Wu, and Resch (1978) on establishing useful a priori prediction algorithms. A 
full report of our work is available in Moran and Penfield (1976). 


BACKGROUND PHYSICS 

The excess propagation path length is given by 


L = 



N(h)dh 


( 1 ) 


where N(h) is the refractivity of the air as a function of height. The refractivity of moist air is de- 
scribed by the Smith-Weintraub equation (Bean and Dutton, 1966), 


N = 


11.6 

T 



4810e\ 
T /’ 


( 2 ) 


where 
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T = temperature (°K) 

P = total pressure (mb) 
e = partial pressure of water vapor (mb). 

The first term in equation (2) arises from the displacement polarizations of all the air constituents 
including about a 1 -percent contribution from water vapor at the surface. It is called the “dry 
term.” The second term in equation (2) is due to the dipole moment of the water vapor molecule 
and is called the “wet term.” 

The refractivity of air is essentially independent of frequency from 0 to 30 GHz. The dispersive 
component of the index of refraction associated with the 22.235 GHz transition of water vapor has 
a refractivity which is less than 0.02 (Liebe, 1969). 

Using the ideal gas law, the dry and wet components of the refractivity, N D and N v , can be written 
as 

N d = = 2.70 x 10 4 p D (3) 

N v = 3.73 x 10 s — = 1720 § , (4) 

v T2 1 

where p D and p v are the densities of dry air and water vapor in grams per cubic meter. At the sur- 
face, N d typically varies between 250 and 300; while in New England, N v varies from between 
about 10 and 100. 

The dry gas obeys the equation of hydrostatic equilibrium which leads immediately to the result 
that the excess dry path is 


Ln = 10 6 


I 


N d dh = AP 0 


( 5 ) 


where P Q is the total pressure at the surface, A = 77.6 — = 0.2276 cm mb -1 , R is the universal 

gas constant, m is the molecular weight of dry air, and g is the surface gravity constant. Hence, the 
value of L d at sea level at the standard pressure of 1013 mb is 231 cm at 45° latitude. L D can 
therefore be estimated to an accuracy of less than 1 cm provided the pressure is measured to an ac- 
curacy of a few millibars since departures from hydrostatic equilibrium are small (Hopfield, 1971). 

The partial pressure and density of water vapor are related by the ideal gas law, so that 


217e _ 3 

gm 3 


Pv = 


T 


( 6 ) 
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The wet path length is therefore 


L v = 1720 



P v ( h ) 
T(h) 


dh. 


( 7 ) 


The assumption that p v is an exponential function with a scale height of 2.2 km and T is constant at 
290° K leads to the approximate formulas 


Ly ~ 

1.3p v 

(8) 

Lv - 

l.Oe. 

(9) 


More exact analysis for prediction of path length from surface meteorological data are given by 
Saastamoinen (1973). 


The brightness temperature at the surface of the earth, derived from the equation of radiative trans- 
fer, can be written as 


T„(iO = T c e + T(£)a(i/,£)e-’-'(^ g )dg, 


where 


and 


T v 



ol{v, £)d£ 


a(v, £)d£ 


( 10 ) 


(ID 

( 12 ) 


and where £ is the distance along the ray path from the observer, a (i>, £) is the absorption coeffi- 
cient, T c is the brightness temperature of any extraterrestial radiation source, r„ is the total atmo- 
spheric opacity along the ray path, and r„ is the opacity between the point of emission and the 
observer. 

The principal contribution to the microwave absorption coefficient are water vapor and oxygen. 
The opacity due to oxygen at 22 GHz is about 0.013 nepers; therefore, the brightness temperature 
contribution is only about 4°K. This contribution changes very little with time. The absorption 
coefficient for oxygen, taken from Meeks and Lilley ( 1 960) has been included in our calculations. 
The absorption coefficient due to the water vapor for the 6 16 -5 23 transition having a rest fre- 
quency, v 0 , of 22.23508 GHz is given by Staelin (1966) as 
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ot(v, £) = 3.24 x 10" 4 e- 664 / T PPv (l +0.0147 

X 3 - 125 ' P 1 


1 


1 


L i v ~ v o ) 2 +A *' 2 C^+^o ) 2 + Ai > 2 J 


+ 2.55 x 10 -8 p y v 2 


Av 

X 3 /2 


cm" 


where 


Av = 2.58 x 10- 3 (l +0.0147 
and v is the frequency in GHz. 


Pv T 


p 1 (T/318) 0 - 625 


The brightness temperature at center line, for the case where t v < 1 , is therefore 


T r <* 


I. 


Pv 

O PT 875 


g-644/T dh 


whereas 


Ly a { 

J n 


Pv 


dh. 


(13) 


(14) 


(15) 


(16) 


To a first approximation, we find that T B (H 2 0, v = 22 GHz) (°K) ~ 2.1 Ly (cm). P decreases by 
10 percent per kilometer and T by 2 percent per kilometer so that a given amount of water vapor 
contributes more heavily at higher altitudes to the brightness temperature than to the wet path 
length although saturation tends to reduce the discrepancy. Frequencies can be chosen to maximize 
the correlation between path length and brightness temperature (Wu, 1979). 


RADIOSONDE STUDIES 

Radiosondes are launched routinely at 1 1 and 23 hours UT from Portland, Maine; Chatham, Massa- 
chusetts; and Albany, New York by the National Weather Service. We obtained data for 50 occa- 
sions during which simultaneous launches were made at the three stations in July and August 1974 
and for 30 occasions in January and February 1974. The data of the significant reporting points 
were used in our analysis. The geopotential height was calculated from the pressure and tempera- 
ture (Hess, 1959), and the water vapor density was calculated from the dew point depression and 
temperature. The absorption coefficient was calculated as a function of height from p v , P, and T. 
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The wet path and the brightness temperatures were calculated from equations (7), (10), (13), and 
(14) using the trapezoidal rule for integration. The mean and rms of Ly were 17±7 cm and 5±2 cm 
for summer and winter. 

The calculated values of wet path length and brightness temperature were compared in various ways 
with a simple least-mean-square analysis program. Some of the results are shown in table 1. To see 
how well the wet path length could be predicted from the surface density of water vapor, p 0 , the 
data were fit to the equation 


Ly = A 1 p Q + A 2 . (17) 

The rms deviations in path length were typically 5 cm for the summer data and 2 cm for the winter 
data. Hence, surface data only offers a slight improvement in the estimate of Ly compared to pre- 
dictions from the climatic mean. 


Table 1 

rms Error in Predicting Wet Path Length from Various Parameters 


Predictor 

summer 

a 

(cm) 

winter 

a 

(cm) 

Seasonal mean ( 1 ) 

6.9 

3.1 

Seasonal mean (2) 

4.4 

2.2 

Surface meteorology (1) 

5.1 

2.1 

Surface meteorology (2) 

3.7 

2.0 

Surface meteorology (3) 

3.9 

2.4 

Remote station radiosonde (1) 

5.0 

1.5 

T 22 and T 19 (1) 

0.3 

0.2 

T 22 (1) 

0.4 

0.4 


NOTES: (1) Portland-Chatham- Albany data from this study. 

(2) Huntington, W. Va. data from Schaper et al. (1970). 

(3) Albany, N.Y. data from Crane (1976). 
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The hypothesis that the wet path at one station might be useful in predicting the wet path at 
another was examined. Wet path data from the station i was fit to the wet path data from the sta- 
tion j with the equation 


L v . = A 1 L v . + A 2 fefcj. (18) 

The mean of the rms deviations were 5 cm for the summer data and 1 .5 cm for the winter. The wet 
path length is therefore fairly well decorrelated over a distance of 200km. Hence, the remote pro- 
files are no better than local measurements of surface meteorology for the prediction of the wet 
path length. 

The correlation between brightness temperatures at 22 and 19 GHz, T 22 and T 19 , and wet path 
length were calculated by fitting the data to the equations 

Ly = AjT 22 (19) 


and 


L v — AjT 22 + A 2 T 19 + A3. (20) 

The rms deviations for the two-temperature fit for the summer and winter data were 0.3 cm and 0.2 
cm, showing that wet path length and brightness temperature are well correlated. 


RADIOMETERS 

The two radiometers in our study were built at the National Radio Astronomy Observatory in 1968. 
Their characteristics are summarized in table 2. In order to measure the antenna temperatures at 
two frequencies, 19 and 22 GHz, each radiometer was configured as a dual-Dicke-switched radi- 
ometer which switched between the antenna and an ambient load at 50 Hz and between frequen- 
cies at 2 Hz. Although the radiometers were very stable, drifting less than 0.3°K (rms) over many 
hours, absolute calibration was difficult. The accuracy of the noise tube calibration was measured 
to be about 1 .5 percent. Hence, the absolute calibration with respect to the 295°K load was about 
3.5°K. To avoid the need for absolute measurements, observations were made at different elevation 
angles. The observable used was the brightness temperature at 30° elevation angle minus that at the 
zenith. For a temperature difference of about 30°K, the calibration error was about 0.5°K. The 
synchronous detectors suffered a peculiar problem that increased the rms accuracy of the tempera- 
ture measurements by 1.5°K. This was, therefore, the dominant contribution to the total measure- 
ment error of about 1 .6°K. 


THE AUGUST 1975 RADIOSONDE EXPERIMENT 

In August 1975, an experiment was performed at the Haystack Observatory in Westford, Massa- 
chusetts, to determine the accuracy of using data from the NRAO radiometers to predict the wet 
path length. A direct comparison of the radiosonde and radiometer data was made at that time. 
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Table 2 

Radiometer Characteristics 


Source 

NRAO (constructed 1968) 

Antenna 

92 cm paraboloid 

Beamwidth (FWHM) 

61 arc minutes 

Peak sidelobe 

-22 db 

Pointing accuracy 

±0.5° (potientometer encoders) 

Receiver 

Superhetrodyne, mixer first stage 

Receiver Temperature 

1500°K 

Bandwidth 

300 MHz 

AT (rms) (4 sec) 

0.1°K 

AT (rms) (5 hours) 

0.2°K 

Noise tube accuracy 

1.5% 

Reference load temperature 

294° K 

Tuning range 

18-24 GHz 

Dicke switching 

Frequency switch @ 2 Hz; load switch @ 50 Hz 

Observation frequencies 

19.0, 22.2 GHz 

Absolute calibration 

4°K 

Relative calibration 
(T(30° elevation) - T(zenith)) 

0.9°K 
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The radiosondes were launched by NWS at the NASA Wallops Island Facility. The radiosondes con- 
tained carbon elements for sensing relative humidity (Brousaides and Morrissey, 1 974). Forty-five 
radiosondes were launched over an 1 1-day period from August 14 to 25, 1975. The basic data were a 
series of measurements of pressure, temperature, and relative humidity from the radiosondes and 
the tracking angles of the telemetry antenna. NWS supplied us with reduced profile data consisting 
of pressure, temperature, water vapor density, wind speed and direction, and other quantities de- 
rived from them, at intervals of 100 feet (after interpolation) from the surface to an altitude of 
30,000 feet. Typically, the radiosondes drifted eastward from Haystack and were at an elevation 
angle of 30° at 1 km altitude. The total pressure, partial pressure of water vapor, and temperature 
at the surface are plotted in figure 1 and show that during the period of the experiment, two weak 
frontal disturbances passed. The mean profile of water vapor density is shown in figure 2. The pro- 
file has, approximately, an exponential dependence with a scale height of 2.2 km, although the high 
altitude part decreases more slowly than an exponential. Individual profiles deviated greatly from 
exponential dependence. 



14 15 16 17 18 19 20 21 22 23 24 25 26 


DAY (AUGUST 1975) 

Figure 1 . The temperature, partial pressure of water vapor, and total 
pressure measured at the surface before each radiosonde launch in 
August 1975. The symbols denote different cloud cover conditions: 
(•) <2/10 cover or thin cirrus (21 launches); (x) >2/10 and <6/10 
cover (13 launches); (+) >6/10 cover (11 launches). 
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Figure 2, The mean profile of the water vapor from the data of 45 
radiosonde launches at Haystack Observatory in August 1975. The 
dotted lines denote plus and minus one standard deviation. 

During each balloon flight, which lasted about 30 minutes, six tipping scans were made at different 
azimuthal angles. The values of brightness temperature, T(30°)-T(zenith), were averaged for each 
balloon flight. The set of 45 brightness temperatures for each frequency and radiometer were com- 
pared to those calculated from the radiosonde data (as described in section II) in order to calibrate 
the radiometers. The temperature data sets (i.e., temperatures at two frequencies from the two 
radiometers and the radiosonde estimates) were fit to the equation 

T i = A 1 T j +A 2 i*i. (21) 

The values of A 2 were typically less than 1°K. The rms deviation was 1.2°K in the comparison of 
data from radiometers 1 and 2, implying the radiometer measurements were accurate to 0.9°K 
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(i.e., 1 ,2°/y/2 ). Most of this error was caused by the problem with synchronous detector. The ex- 
pected radiometer error was 1 .6/^/6 or 0.7°K. The rms deviation in the fit of radiometer temper- 
atures to radiosonde predicted temperatures was typically 2.5°K. This suggests that most of the 
experimental error was in the radiosonde data, or due to the effects of clouds. 

Extensive correlation analysis was performed between radiometer temperatures and the wet path 
and theoretical temperatures estimated from radiosonde data. A plot of the wet path versus bright- 
ness temperature at 22 GHz is shown in figure 3. The path length predicted by radiometer 1 from 
equation (20) is plotted versus the path length estimated from the radiosonde data in figure 4. 



Figure 3. The wet path lengths for the vertical direction 
calculated from radiosonde profiles versus the brightness 
temperature (T(30°)-T (zenith)) at 22.2 GHz measured by 
one of the NRAO radiometers. The symbols, denoting dif- 
ferent cloud conditions, are defined in the caption to figure 
1. The non-zero intercept is due largely to oxygen emission. 


CONCLUSIONS 

The residuals from various regression analyses for the August 1975 experiment are listed in table 3. 
We reached the following conclusions: 

1 . The mean value of the path lengths estimated from radiosonde data for the 45 launches 
was 14.8 cm with an rms deviation of 5.6 cm. For the 21 radiosondes launched in clear or 
nearly clear skies, the mean was 10.3 ± 2.5 cm. 
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Figure 4. The wet path length, estimated from the formula 
Ly = T 1 T 22 + A 2 T 19 + A 3 with temperatures from radiom- 
eter 1 versus the wet path lengths derived from radiosonde 
data. The coefficients A 1 , A 2 , and A 3 were derived by 
fitting the brightness temperatures from the radiometers to 
the wet path length estimated from the radiosonde data. 

2. Prediction of the wet path length from the surface partial pressure of H 2 0 had an rms 
deviation of 3.2 cm, an improvement over the rms deviation about the mean. 

3. Prediction of the wet path length at Haystack from those estimated for Portland, Albany, 
and Chatham gave no better results than prediction from local surface meteorology. 

4. The rms deviation in wet path length estimated from radiometer data was 1 .5 cm. 

5. The rms deviation in path length predicted from brightness temperatures computed from 
radiosonde data was only 0.4 cm. This may represent the ultimate accuracy of the tech- 
nique without adding higher-order regression coefficients or correcting for the effects sat- 
uration caused by large microwave opacity. 


The results from this experiment were consistent with the results from our previous analysis of 
radiosondes launched at Albany, Portland, and Chatham. The residuals for the winter data were, 
of course, much smaller. 
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Table 3 

Results of Regression Analysis on August 1975 Data 


RMS Residuals in Wet Path Length 


Predictor 

FitG> 

All dataH) 

Best Weather data! 2 ) 

Mean 


5.6 cm 

2.5 cm 

Surface partial pressure 

Ly Aj C "t“ A 2 

3.2 

1.5 

Path length from radio- 
sondes at Portland, 
Albany, Chatham 

Ly - Al Lyj + A 2 Ly 2 + A 3 Ly 3 

3.0 


Brightness temper- 
atures from NRAO 
radiometers 

Lv = AiT 2 2 “1 A . 2 Tu + A 3 

1.5 

1.2 

Brightness temper- 
atures estimated from 
radiosonde data 

L v = Ai T 22 + A 2 T 19 + A 3 

0.4 

0.2 


NOTES: (1) 45 data points. 

(2) 22 data points for which cloud cover was <0.2. 

(3) Ly estimated from radiosonde data. 


The stability of the regression coefficients derived presents a problem requiring some study. Wu 
(1979) has shown that coefficients can be found which are insensitive to climatic conditions if cor- 
rect pairs of frequencies are chosen. Wu requires that the coefficients satisfy the relation (A l /A 2 ) = 
~(v 2 /v 1 ) 2 so that the brightness temperature of clouds, which is proportional to v 2 , does not effect 
the path length estimate. With this constraint, 19 and 22 GHz are not a good frequency pair. A 
good way to include the radiometer data in the analysis of VLBI data may be to add the observed 
brightness temperatures to the data base of measured interferometer observables and estimate the 
regression coefficients for each site in the general data analysis. In this way the stability of the coef- 
ficients under different climatic conditions is less important. 

This work was supported under NASA contract NAS5-20975. We thank Hays Penfield and Joe 
Hayes for construction and operating support, and 1. 1. Shapiro, C. C. Counselman, and C. A. Knight 
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